Inverse Spin Hall Effect by Spin Injection 
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Motivated by a recent experiment [Nature 442, 176 (2006)], we present a quantitative microscopic 
theory to investigate the inverse spin-Hall effect with spin injection into aluminum considering 
both intrinsic and extrinsic spin-orbit couplings using the orthogonalized-plane-wave method. Our 
theoretical results are in good agreement with the experimental data. It is also clear that the 
magnitude of the anomalous Hall resistivity is mainly due to contributions from extrinsic skew 
scattering, while its spatial variation is determined by the intrinsic spin-orbit coupling. 



PACS numbers: 72.25.Ba, 72.25.Hg, 71.70.Ej 



In the presence of an applied electric field, spin-orbit 
(SO) coupling can cause deflection of electrons with op- 
posite spins in opposite directions, resulting in a nonvan- 
ishing transverse spin current. This so-called spin-Hall 
effect (SHE), which was first predicted by D'yakonov 
and Pcrel in 1971 pj and was recently emphasized by 
Hirsch[|, has attracted a great deal of experimental and 
theoretical interest Q because of its potential applica- 
tions in Spintronics. Depending upon its origin, the SHE 
is classified into two types: intrinsic [l], [H, 3 and extrinsic 

she nn. 

Experimental observations of SHE have also 
been reported recently 0, @, E3] • 

The SO coupling mechanism can also lead to a recip- 
rocal SHE effect, the "inverse spin-Hall effect" (ISHE), 
which was predicted by HirschQ ■ When a purely longitu- 
dinal spin current is applied, the electrons with opposite 
spins, flowing in opposing longitudinal directions, move 
toward the same transverse side of the sample due to SO 
interaction, resulting in charge accumulation. Recently, 
Valenzuela and Tinkham presented a first clear ISHE 
observation 11|. By a nonlocal spin injection into metal- 
lic aluminum through a ferromagnetic/nonmagnetic con- 
tact (the configuration is illustrated in Fig. 1), they ob- 
served a finite Hall resistivity of order of milli-Ohms over 
a spin-diffusion length of the order of /im. 

It is well known that spin[I3] and charge [r| relax- 
ations in aluminum are strongly influenced by band- 
structure anomalies which arise mainly from several acci- 
dental degeneracy points (near W points) in the Brillouin 
zone (BZ). Therefore, in an analysis of the ISHE by spin 
injection in aluminum, the momentum dependence of the 
scattering should be taken into account. Earlier work by 
Zhang analyzing the SHE in the presence of spin diffu- 
sion treated only momentum-independent relaxation [lij. 
In this Letter, we provides a realistic, quantitative the- 
ory for ISHE by spin injection in aluminum consider- 
ing the momentum dependencies of the scattering rates 
and of the side-jump and skew scattering contributions 
to the anomalous Hall current (AHC). Both the intrin- 
sic and extrinsic SO couplings are taken into account 
and the electronic bands and states are determined by 
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FIG. 1: (Color online) Measurement apparatus configuration 
for the spin injection experiment of Ref. A net current 
/total is injected through aluminum (Al) via a ferromagnetic 
contact and it flows away from the Hall measurement region. 
Near the contact, i.e. at x = 0, /total = + Ili, Ili = 
Il(x — 0), and If(x = 0) + ij.(a; = 0) = 0, resulting in 
a vanishing net electric current for x > and a finite spin 
current, If(x) — Ii(x), over a distance of order of the spin- 
diffusion length. 



the orthogonalized-plane-wave (OPW) method. The re- 
sulting Hall resistivities are in good agreement with the 
experimental data. It is also clear that, with spin in- 
jection into aluminum, the spatial variation of AHC is 
determined by the intrinsic SO coupling while its mag- 
nitude is mainly due to contributions from the extrinsic 
skew scattering. 

In general, the SO interaction induced by ionic fields, 
i.e. the intrinsic SO coupling, may affect both the en- 
ergy bands and the states of electrons in solids. How- 
ever, in aluminum, its effect on the electronic band 
structure can be ignored. Hence, in our treatment, 
the energy band structure of electrons is determined by 
the standard OPW method in the absence of SO in- 
teraction, while the electronic Bloch states are a mix- 
ture of spin up and spin down species: ipfip (r) = 



a$(G)\ t i>+b%>(G)\}i> 

fi = 1, 2 representing spinors (p, = 3— /x), p as lattice mo- 
mentum confined to the first BZ, n as band index, and G 
as reciprocal lattice vector. The coefficient a„p^(G) is de- 



exp[z(p — G) • r], with 



2 



termined by the pseudopotential form of the Schrodinger 
equation, while fr^p (G) is obtained from first-order per- 
turbation theory by taking account of the SO part of the 
pseudopotential, AL • SVi (A is the SO coupling constant 
and V\ is the operator projecting on the orbital angular 
momentum state 1 = 1). 

In the spirit of the pseudopotential method, the 
electron-impurity scattering potential, V(p, k), corre- 
sponding to scattering of an electron from state p to state 
k, is the difference between the impurity- and host-atom 
form factors with a structure factor correction for lattice 
distortion around the substitutional impurities. It in- 
volves not only a spin- independent part Vb(p, k) but also 
a SO part y so (p,k): V(p,k) = V (p, k) + V so (p, k). In 
general, the elements of Vso(p,k), Vg^(p, k), take the 
form, Vg^(p, k) = A(p, k)p x k- < v\<r\(i > with a as 
the Pauli matrices, and A(p, k) depending only on the 
magnitudes of the electron momenta p and k. 

In the presence of spin diffusion, the transport proper- 
ties of solids may be described by a Wigner distribution 
function p(n, p, R), which is a 2 x 2 diagonal matrix de- 
pending on the band index n, and on both Wigner coor- 
dinates: microscopic electron momentum p and macro- 
scopic space coordinate R = (x,y,z). Under steady state 
conditions, p(n, p, R) obeys a kinetic (Boltzmann) equa- 
tion given by 

£.V R p wt (n,p,R)= E {^T'(P-k) 

p,k,n' 



x \Pw(n', k, R) - Pw(n, p, R)' 
+W$'(p, k)^ p (n', k, R) - p w (n, p, R)]} 



(1) 



Here, W™" (p, k) and W™^ (p, k), respectively, are the 
non-spin-flip and spin-flip scattering rates. They can be 
obtained by analyzing the matrix elements of the scat- 
tering potential: 



W™ (p,k) = 2ir6(e np - e„, k ) 



E {K^(G') + ^L(G') 

[< ! (G) + ^(G)] V (p- G,k- G') 
[all(GOraW(G)V^(p- G,k- G')}f (2) 



G.G 

X 



and 



W™ (p,k) = 2n6(e np ~ s n , k ) 



E {(< } (G)[6^(G')]* 

+^ } (G)[a^(G')]*) V (p - G,k- G') 
+a^(G)[oi« (G')]*V^(p G,k G')}| 2 , (3) 



with e np as the energy dispersion relation for n-th band 
electrons. Note that in Eqs. (2) and (3) both the in- 
trinsic and extrinsic SO couplings are involved. The 
last term on the right-hand side of Eq. ([T]) leads to an 
exchange of electrons with different spins, causing spin 
relaxation. This so-called Elliot- Yafet spin-relaxation 
mechanism (lit flif in aluminum has been carefully in- 
vestigated by Fabian and Das Sarma taking account of 
electron-phonon scattering [12]. They found that spin re- 
laxation in aluminum is determined mainly by several 
accidental degeneracy points near the W points in BZ. 



It has already been demonstrated that intrinsic SO 
coupling can make a nonvanishing contribution to AHC if 
the off-diagonal elements of the distribution function are 
finite 17j. However, in the spin injection experiment, the 
distribution function is essentially diagonal since the off- 
diagonal driving forces, such as the electric field and Zee- 
man energy splitting, are negligible. In particular, in the 
experiment of Ref. , there is no external driving elec- 
tric field in the Hall measurement region, while the Zee- 
man energy splitting for the maximum applied magnetic 
field, B = 3.5T, is only 0.2meV, much smaller than the 
free-electron Fermi energy of aluminum, = 11.7 eV. 
Hence, to determine AHC, one only needs to consider the 
extrinsic SO coupling. 



The extrinsic SO interaction makes nonvanishing con- 
tributions to AHC through two mechanisms: a side-jump 
process proposed by Bergerfljj] and also skew scattering 
given by Smit [l^| . The side-jump AHC arises from a 
sidewise shift of the center of the electron wave packet, 
while the skew scattering contribution corresponds to 
an anisotropic enhancement of the wave packet due to 
electron-impurity scattering described in the second Born 
approximation |20j. Obviously, the analysis of spin injec- 
tion in aluminum calls for the formulation of these two 
mechanisms of AHC in terms of a distribution function. 
In previous studies, these two AHC mechanisms were de- 
scribed in the framework of a momentum-independent 
relaxation time[2ll. [22l|. or in linear response theory(23j|. 
and hence the existing results can not be applied in the 
present study of the ISHE by spin injection in aluminum. 



To analyze the side-jump process, we begin with the 
anomalous term of the current operator, which is as- 
sociated with function < ip+(n, i)^ M (r, t) > (Vv(r, t) 
and tjj^ (r, t) are electron field annihilation and creation 
operators, respectively). Expanding the statistical av- 
erage involved in this function to first order of the 
electron-impurity interaction and applying the Langreth 
algebra 24j taken jointly with the generalized Kadanoff- 
Baym ansatz 2|| 26 1, we finally arrive at the form of the 
steady-state current, J S,I (R), in the lowest-order gradi- 
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ent expansion of the Wigner coordinate R: 
Jf (R) = -nreN t £ £ (-l)^W(np, n'k, G x , G 2 ) 



E 

p,q Gi,g 2 

ju,n,n G 3 ,G 4 



x? W (n'k, np, G 3 , G 4 )e a0z [(k - G 2 ) - (p - d)]^ 
x<5(e„'k - e„ p )[/3 /iAl (n,p,R) -/5 w (n',k, R)], 



with 



(GOJ'o^CGOAdp - GxlJk G 2 



A^)(np,n'k,Gi,G 2 ) 



F W (np,n'k,Gi,G 2 ) = 



(G 2 )]*a$(Gi)V (T> - 



Gi,k — G 2 ), and s a p z as the totally antisymmetric 
tensor. 

The skew scattering AHC, J SS (R), is associated with 
the additional SO term of the distribution function in the 
second Born approximation. Substituting this term into 
the non-SO part of current operator, in the first order of 
SO coupling constant A, J SS (R) takes the form, 



the 15 shortest reciprocal lattice vectors are considered 
but the pseudopotentials, except for Ashcroft's Vm and 
^200 113 1 are assumed to vanish. This treatment enables 
us to obtain a continuous electron velocity when the elec- 
tron momentum crosses from one symmetry section of 
the fee BZ into another. The sums over all electron 
(4) momenta in the first BZ are performed by the tetrahe- 
dron method[28]. Further, we assume that the electron- 
impurity interaction arises from defects, which may be 
introduced in the process of material fabrication. In 
this way, the scattering potential is just given by the 
form factor of aluminum but with the opposite sign. In 
the calculations, the spin-independent part of scatter- 
ing potential is described by a screened Ashcroft form 
factor [29], while the SO part is chosen from Ref.[I(| with 



J SS (R) = i^eNi 
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:P/x/x(n',k,R)r(n',k), 



(5) 



with v„ p = d£ np /dp as the electron velocity and Ni as 
the impurity density. r(n',k) is the transport relaxation 
time, which has been carefully investigated in Ref. [l3T |. 
The last two terms in curly brackets of Eq. ([5]) are ob- 
tained from the first term by following the replacement 
rules indicated in the square brackets. 

It should be noted that the side-jump and skew scat- 
tering contributions to AHC, J„ J (R) and J^ S (R), ex- 
ist in conjunction with a nonvanishing longitudinal spin 
current, rather than the net electric current. In a sys- 
tem with a longitudinal charge current but without mag- 
netization, the electrons with opposite spins distribute 
equally and the factor (— 1) M in Eqs. Q and (O leads 
to vanishing of the AHC. However, when a purely longi- 
tudinal spin current is present, we have pj^ (n', k, R) = 
— y3nn (n', k, R) with pW(n',k,R) as the deviation from 
the equilibrium Fermi distribution. Hence, the signs of 
contributions to Hall current from electrons having op- 
posite spins are the same, resulting in a finite transverse 
charge current. 

We perform a numerical calculation to investigate 
the ISHE in the spin injection experiment illustrated 
in Fig. 1. The electronic band structure of aluminum 
is evaluated by a modified 4 OPW method, in which 



SO coupling constant A = 2.4 x 10 
tcrmine the transverse Hall current, 



d a.u.[ia. To de- 

J' a 



tSJ 

J v 



tSS 

•Jy , 



and the longitudinal spin current, J%(x) = J] x {x) 
Jlx(x) with J^ x (x) = Y^,n,p e ( v np)xP^(p,x), we solve 
the Boltzmann equation, Eq. fl}, using the Fermi- 
surface-harmonic-expansion method (3l|. assuming an 
initial deviation from the equilibrium Fermi distribu- 
tion, pfy(n,p,x = 0), of form: pfy(n,p,x — 0) = 
(— l) M eE c • v„ p T np np(£„p) with n' F (e) as the derivative 
of the Fermi function. The effective electric field E c is 
obtained by considering the initial spin-up current densi- 
ties injected into the right-hand side of the aluminum at 
x = in the experiment of Ref. 1JJ: J-[ x {x = 0) = 3.75 
or 1.8 x 10 9 A/m 2 for samples with thicknesses t = 12 
or 25 nm (in the experiment, /total = 50 fiA, the sam- 
ple widths are w — 400 nm and the spin-polarization, 
P = [I L ] - /Li)//totai at x = 0, is 0.28). The impu- 
rity densities are determined from the mobilities of the 
experimental samples. 

In Fig. 2, we plot the longitudinal spin current 
and the AHC J y as functions of coordinate x for the 
two different values of J^ x (x = 0). It is evident that 
with increasing x, J x (x) and J y (x) decrease exponen- 
tially: J x (x), Jy(x) ~ exp(— x/x s d) with x s d as the spin- 
diffusion length. For aluminum, we find ^'^ oory = 830 nm 
for t = 25 nm, and 2^ c ° ry = 512 nm for t = 12 nm, which 
are in good agreement with the experimental results: the 
experimental x s d values are = 735 nm for t = 25 nm 
and = 490 nm for t = 12 nm. 

From Fig. 2 we see that the AHC exists but is rela- 
tively small: it is almost five orders of magnitude smaller 
than J x (x). Nevertheless, the corresponding Hall volt- 
age can be detected. We calculated the Hall resistivities 
for the two initial values of J-\ x {x) and compared them 
with the experimental data. The results are plotted in 
Fig. 3. Our theoretical results, without any adjustable 
parameters, are in good agreement with the data. For 
t = 12nm (or Jf K (0) = 3750A/m 2 ), our theoretical 
values for -Rahe agree perfectly with the data. For 
t = 25 nm (or J Tx (0) = 1800 A/m 2 ), the calculated re- 
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FIG. 2: (Color online) Coordinate dependencies of the longi- 
tudinal spin current, J* (a), and anomalous Hall current, J y 
(b), in the spin injection experiment for different initial values 
of J u {x): J u (0) = 3750, and 1800 A/^m 2 . 

suits are a bit larger than those obtained experimentally. 
This may be due to the anisotropy of the experimental 
sample, which could lead to a different A-value than the 
one we use here. Note that the sign of the calculated Hall 
voltage also agrees with that found experimentally: the 
positive spin current, J^(x), polarized along the positive 
^-direction produces a negative Hall electric field. 

Our numerical calculation also shows that the contri- 
butions to AHC are dominated by skew scattering. The 
side-jump AHC contribution is only about 3.5% of the 
skew scattering one. This is due to the fact that the mo- 
bilities of the samples involved are small and the impuri- 
ties are relatively dense: Ni is about 1% of the electron 
density. Note that the signs of both the side-jump and 
skew scattering AHC contributions are the same. 

We also find that the effect of the extrinsic SO coupling 
on the diffusion length vanishes. This holds valid for 
materials with space inversion symmetry (SIS). In the 
case of weak SO coupling, a^p is real: atfp = [a«p]*, and 
SIS leads to alp' = a^_ p = alp'. Hence, the extrinsic 
parts of the scattering rates are asymmetric under the 
momentum inversion operation, resulting in a vanishing 
contribution to the diffusion length. This implies that 
the spatial dependence of AHC is determined only by the 
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FIG. 3: (Color online) Anomalous Hall resistivities, 7?ahe, as 
functions of coordinate x for different initial values of J^ x (x): 
Jtz(0) = 3750, 1800 A//im 2 . The solid and dotted lines are 
our theoretical results, and the squares (circles) are the ex- 
perimental data for Al thicknesses t = 12 (25) nm (Ref. |ll|). 



intrinsic SO coupling since J^ J (R) and Jf s (R) depend 
on R only through the distribution function. 

In conclusion, we have presented a fully microscopic 
theory to analyze the ISHE observed in the spin-injection 
experiment in aluminum taking account of both the in- 
trinsic and extrinsic SO couplings. Employing the OPW 
method, the momentum-dependent scattering rates of 
the Boltzmann equation were treated carefully and the 
side-jump and skew scattering contributions to the AHC 
were formulated in terms of the distribution function. 
Performing realistic calculations for aluminum, we have 
determined anomalous Hall resistivities that are in good 
agreement with the experimental data. 
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